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1. Introduction 



The potential formed by two heavy quarks is among the early applications of Quantum Chro- 
modynamics (QCD). It is an important ingredient in the description of the properties of heavy- 
quark bound states. Among the most prominent applications, which require high-order results of 
the static potential, are the top quark production cross section close to threshold and the extraction 
of the bottom quark mass from T sum rules (see, e.g., Ref. ^ for a review). 

The n-loop corrections to the quark anti-quark potential are usually parameterized by the con- 
stants a, such that in momentum space it takes the form 
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where the renormalization scale has been identified with \q\ and for SU{Nc) we have Ca 
Cf = (A^2 _ \)/{2Nc) and A^^ = 3. 

The one-loop corrections to have been computed more than 30 years ago [Q]. Around 

the same time it has been discovered that V(|^|) is not infra-red safe and starting from three-loop 
order divergences appears [||]. In the 1990ies the two-loop corrections have been computed in the 
works |5|] and shortly afterwards the coefficient of the three-loop divergence has been evaluated 
indirectly by examining the ultra-soft contribution to the energy of two heavy quarks ^ 0]. First 
steps towards the finite part of the three-loop corrections have been performed in Refs. [||, ^ where 
the fermionic contributions have been evaluated. The gluonic contribution which completes the 



knowledge about has been obtained by two independent computations [ ]10| , 

The evaluation of Refs. ||^ |l^ is highly automated in order to avoid errors due to manual 



interactions. After generating the amplitudes for the Feynman diagrams with QGRAF [ ]12| ] we use 
q2e and exp |l^] in order to rewrite the expressions in FORM [15] format which is used for 
taking the traces and further simplifications. The reduction to master integrals is performed with 



the program package FIRE [16] and the resulting master integrals are computed with the help of 



the Mellin-Barnes technique (see, e.g., Refs. [17, 18, 19|). Cross checks are based on the program 



FIESTA [gOp which incorporates the sector decomposition algorithm. 

We managed to compute all the necessary coefficients of the e expansion of the master inte- 
grals analytically with the exception of three terms of order e ^ . Results for some master integrals 
are presented in the next section. Let us close this Section by summarizing the results for ai, from 
Refs. [§, 10]: It is convenient to decompose as according to the powers of «/, the numbers of 
massless quarks = af^^ + af^n^ + a'^^'ni + af^ . These are our results: 
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jabcd jahcd 

af = 502.24(1) C] - 136.39(12) ^ , (1.2) 

where Tp = 1 /2, df^^^df'^ /Na = (18 - 6N^+N^) /{96N^) and df'^^df"^ /Na = {N^ + 6A^c) /48. 



2. Results for selected master integrals 



If one counts all master integrals for all the types of the integrals appearing in the calculation 
in the general i§ -gauge one obtains around hundred master integrals. However, only 41 master 
integrals contribute to the three-loop static potential in the Feynman gauge. Eighteen most compli- 
cated master integrals are shown in Fig. |l|. Only two of them, /n and /is, are non-planar. Diagram 
/i7 does not involve static lines and is just a three-loop propagator diagram which can be found in 



Ref. 1 21] and which we need up to order e^. Diagram /ig represents the most complicated master 
integral in our calculation. Explicidy, we have /ig 



^^W,o,i,0' where 
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with the causal —iO implied in all propagators. It turns out that it is more convenient to use, instead 
of this master integral, a master integral with a numerator: /ig = F^^^l -2 i o- ^^^^ finite at d = 4 
and its value is one of three (yet) analytically unknown constants. The other two pieces of the 
three-loop static potential which are only known numerically are the 0{e) terms in In and /ig. 




Figure 1: Most complicated master integrals contributing to a^. The solid and zig-zag lines correspond to 
relativistic and static propagators, respectively. Small black boxes denote monomials in numerators. 
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In the following we present our analytic results for the 16 most compUcated master integrals 
which are conveniently expressed as special cases of one of the four functions 
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where in all propagators, apart #9, the causal —iO is impUed. Here are our results: 
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where Q = \J —(f- and (jTT^'/^e"^^)^ is implied as a factor on the right-hand side. These e- 
expansions are up to the order which contributes to the static potential, with the exception of l\ \ 
and /16, where one more order is desirable. Results for all the master integrals as well as details of 
their calculation will be published elsewhere. 
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